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We present the first results of numerical test of the Kugo-Ojima confinement criterion in the lattice Landau 
gauge. The Kugo-Ojima criterion of colour confinement in the BRS formulation of the continuum gauge theory 
is given by ul{0) — —S^, where 

die''"^-^'(0|rD^c"(x)3(A, X cfiym = (g,, - P^)utip^). (*) 

We measured the lattice version of 11^(0) in use of l/{—dD{A)) where D^{A) is a lattice covariant derivative in 
the new definition of the gauge field as f/ = e^. We obtained that ul{0) is consistent with —c5^, c = 0.7 in SU(3) 
quenched simulation data oi (3 — 5.5, on 8* and 12^. We report the /3 dependence and finite-size effect of c. 



1. INTRODUCTION 

The colour confinement problem in the contin- 
uum gauge theory was extensively analysed in use 
of the BRS formulation by Kugo and Ojima[Q. 

The QCD lagrangian is invariant under the 
BRS transformation and the physical space is 
specified as the one that satisfies the condition 
Vphys = {\phys)} 

QB\phys) = 0. 



d^F^^ + {Q-Bj-D^c}, where its ambiguity by di- 
vergence of antisymmetric tensor should be un- 
derstood, and this ambiguity is utilised so that 
massless contribution may be eliminated for the 
charge, to be well defined. 

Denoting g{A^ x c)° u'^df^j'^, and then 
(1 + u)^9p7'', one obtains the eq.(*) pro- 
has a vanishing expectation value. The 
current {QsjD^c} contains the massless com- 
ponent, (1 + u)^a^/3^(a;). We can modify the 
Noether current for colour charge such that 
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gJ'u = gJ, 



{Qb.D^c). 



and {F X GY = fabcF'G^. 

Under the assumption that BRS singlets 
have positive metric, it is proved that Vphys 
has positive semidefinite in such a way that BRS 
quartet particles appear only in zero norm. 

One finds from the BRS transformation that 
for each colour a, a set of massless asymptotic 
fields x", 7°, 7 form a BRS quartet. 

The Noether current corresponding to the con- 
servation of the colour symmetry is gJ'^ = 



In the case of 1 + u = 0, massless component in 
g Jq is vanishing and the colour charge 



d^'xiQB.g-^Doc^ix)} 



(1) 



becomes well defined. 

The physical state condition Q b \phys) ~ to- 
gether with the equation implies that all BRS 
singlet one particle states |/) G Vphys are colour 
singlet states. This statement implies that all 
coloured particles in Vphys belong to BRS quartet 
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and have zero norm. This is the colour confine- 
ment. 

2. LATTICE CALCULATION OF 

The Faddeev-Popov operator is 

M[U]^-{d-D{A)) = -{D{A)-d), (2) 

where the new definition of the gauge field is 
adopted as U = e"^, and the lattice covariant 
derivative D^{A) = + Ad{Afj) is given in l^. 

The inverse, M^^[U] = (Afo-Mi[J7])-i, is cal- 
culated perturbatively by using the Green func- 
tion of the Poisson equation M^^ = (—9^)^^ and 
Ml = d^,Ad{A^(x)), as 

= + ^ {M-^MifM^^. (3) 

fe=0 

In use of colour source \\°'x) normalised as 
Tr (A*^ a; I A''a;o} = 5°-^5x,xa^ the ghost propagator 
is given by 

G'^\x,y) = {Tr{\''x\{M[U]r^\\'y)) (4) 

where the outmost () specifies average over sam- 
ples U. 

The ghost propagators of /3 = 5 and 5.5 are 
almost the same and they are infrared divergent 

which can be parameterised as . We observed 

that the ghost propagators of /3 = 6 is similar to 
that of /3 = 5.5 and its finite-size effect is small |3|. 

In the similar way, one can calculate the Kugo- 
Ojima parameter at p = as, 

= {Tr{Xy\D^{A){M[U])-\Ad{AM\^p))\,=^ (5) 

We observed that off-diagonal element of is 
consistent to zero, but there are statistical fluc- 
tuations. The projection operator g^i, — ^'^f " in 
equation (*) is treated such that it has an expec- 
tation value I in the limit of 0. 

Making the accuracy of the covariant Laplacian 
equation solver higher, we observe the tendency 
that the expectation value of \u^\ increases. 

At /3 = 8, direct measurement of gives a 
large fluctuation, but suitable twisting treat- 
ment for each sample so that the Polyakov scat- 
ter plot should be concentrated around argz — 0, 
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Figure 1. The dependence of space and 
colour diagonal part of the Kugo-Ojima pa- 
rameter it° on /3 and lattice size. (An aver- 
age over the four directions and eight adjoint 
representations.) The data points are (3 = 
5.5, 8^ X 16;/3 = 5.5,124;/? ^ 5.5, 8^;/? = 
6,124;/3 = 6,84(no Za); P = 6,8^{with Z3); (3 = 
6,8'^{with Z^, minimum Landau); (3 = 
8, 8* {with Z3); /3 = 8, 8'^ {no Z3) respectively from 
left to right. 



suppresses the fluctuation and makes the quality 
of the data better. We consider that this treat- 
ment is indispensable in the simulation where Z3 
symmetry persists and the Z3 factor affects the 
observed quantity. The similar behaviour is ob- 
served in /3 = 6, 8^ lattice. The minimum Landau 
gauge fixing via smeared gauge fixing performed 
at /3 = 6, 8"* lattice does not change the expec- 
tation value obtained after the Z3 twisting but 
reduces the standard deviation. 

The absolute value of u° is plotted as the func- 
tion of the spatial extent of the lattice aL where 
a is calculated by assuming = lOOMeV. We 
find for aL < 2/m, there exists large finite-size 
effectDWe expect that by making L large and a 
small, such that aL > 2/m, the absolute value of 
M° becomes closer to 1. 

Non-symmetric lattice 8^^ x 16 yields non- 
symmetric data in /i of (*). This fact shows neces- 
sity of tuning lattice constants according to the 
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Table 1 



Kugo-Ojima parameter Space-diagonal(/i = i/) and ofF-diagonal components. All data of 8^ are the 
average of 100 samples. 'Z3' and 'min' means Z3 twisting and the minimmn Landau gauge fixing. 









diag 


off - diag 


diagi 


diag2 


diags 


diag 4 


p = 


5.5, 


8^ X 16 


-0.739(135) 


0.002(60) 


-0.776(109) 


-0.779(105) 


-0.818(118) 


-0.581(49) 


f3 = 


5.5, 


12* 


-0.715(46) 


0.003(32) 


-0.729(60) 


-0.713(43) 


-0.705(39) 


-0.712(38) 


(3 = 


5.5, 


84 


-0.664(69) 


0.002(45) 


-0.669(71) 


-0.656(70) 


-0.667(67) 


-0.664(67) 


/3 = 


6.0, 12^ 


-0.548(133) 


-0.015(85) 


-0.555(123) 


-0.561(107) 


-0.508(133) 


-0.566(159) 


/3 = 


6.0, 


8 , with Z'i 


-0.303(80) 


0.002(29) 


-0.286(76) 


-0.307(66) 


-0.325(81) 


-0.293(91) 


P = 


6.0, 


8^, with Z3, min 


-0.308(88) 


-0.000(35) 


-0.312(123) 


-0.311(78) 


-0.317(75) 


-0.292(59) 


P = 


6.0, 


8^, no Z3 


-0.354(176) 


-0.001(76) 


-0.339(130) 


-0.347(161) 


-0.378(239) 


-0.353(151) 


/9- 


8.0, 


8^, with Z3 


-0.183(74) 


0.002(20) 


-0.177(71) 


-0.197(77) 


-0.221(83) 


-0.138(19) 


/3 = 


8.0, 


8'*, no Z3 


-0.338(513) 


0.0116(251) 


-0.264(278) 


-0.359(553) 


-0.334(610) 


-0.394(536) 



Figure 2. The finite-size effect of the Kugo-Ojima 
parameter \u'^\ as the function of the spatial ex- 
tent of the lattice aL{fm). 



non-symmetric lattice size and the lattice dynam- 
ics. 



value decreases as increases. 

We observed the gluon propagator is infrared 
finite H] and the ghost propagator is infrared di- 
vergent, suggested to be more singular than — , 
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but less singular than — r. These results quali- 
P . 

tatively agree with the Gribov-Zwanziger's con- 
jecture|^,||, and are consistent with the results 
of Dyson- Schwinger equation Q. It is nice to 
observe that the infrared finiteness of the gluon 
propagator is in accordance with the Kugo-Ojima 
colour confinement. As stated in their inverse 
Higgs mechanism theorem, if we have no mass- 
less vector poles in all channels of the gauge field, 
AIj^, and if the colour symmetry is not broken at 
aU, it follows that 1 -t- u = 0.0. 

This work is supported by High Energy Accel- 
erator Research Organization, KEK Supercom- 
puter Project(Project No. 99-46), and by Japan 
Society for the Promotion of Science, Grant-in- 
aid for Scientific Research(C) (No.ll640251). 
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